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Abstract 

We have analyzed the D — > ttK, KK decay with the naive factorization 
(NF), QCD factorization (QCDF)and QCD factorization including soft-gluon 
exchanges (QCDF+SGE). In these decay channels, the soft-gluon effects are 
firstly calculated with light cone QCD sum rules. Comparing the three kind 
approaches, we can find the calculation results have made much more im- 
proved QCD factorization (QCDF) than the naive factorization (NF), and the 
calculation results have also made improved QCD factorization including soft- 
gluon exchanges (QCDF+SGE) than the QCD factorization (QCDF) in the 
color-suppressed decay channels. In addition, we find the soft-gluon effects 
are larger than the leading order contributions, and the calculation results 
are close to the experimental data for the color-suppressed decay channels. 
In color-allowed decay channel D° — > tt + K~, the soft-gluon effects are small 
and we should consider other power terms, such as final state interaction and 
annihilation effects. 
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1. Introduction 



The study of heavy meson decays is important for understanding the standard model 
(SM) and search for the sources of CP violation. However, The hadronic two-body weak 
decays of D meson involve nonperturbative strong interactions and spoil the simplicity 
of the short distance behavior of weak interactions. Therefore, a simplified approach 
in which the amplitudes of these processes are given by a factorizable short distance 
current-current effective Hamiltonian is not expected to work well. Various approaches 
were employed to include long distance effects. The most commonly and very frequently 
used prescription, motivated by jj- arguments [1], is to apply generalized factorization 
[2-3]. This phenomenological treatment works reasonably well in color-allowed D decays 
[3], but it is failing in the color-suppressed D — > ixix, nK and D — >■ KK decays [4]. 

It is necessary that we study D meson nonleptonic decays beyond the factorization 
approach. A few years ago, M. Beneke et al. [5] gave a NLO calculation the hadronic matrix 
element of B — > nir, nK in the heavy quark limit. They pointed out that in the heavy 
quark limit the radiative corrections at the order of a s can be calculated with perturbative 
QCD method. In D — > irK decay, the momentum transition square is q 2 = 1.7GeV 2 , 
and the radiative corrections of the hard-gluon exchanges can also be calculated with 
perturbative QCD approach. So, the hadronic matrix elements for D — >■ nK can be 
expanded by the powers of a s and A ^ CD as follows: 

(Ki r \O i \D) = {K\j 1 \D){n\j 2 \0)[l + ^r n a: + O(^)i (1) 

ui c 

where 0{ are some local four-quark operators in the weak effective Hamiltonian and 
are bilinear quark currents. In Eq. (1), the power correct term O ( C ^ D ) includes soft- 
gluon effects, final state interaction, which can not be calculated in QCD factorization 
and perturbative QCD method. For the B meson two-body decay, the term is small, 
but it is large and can not be neglected in the D — > nK decay. A few years ago, A. 
Khodjamirian [6] has presented a new method to calculate the hadronic matrix elements of 
nonleptonic B meson decays within the framework of the light cone QCD sum rules, where 
the nonfactorizable soft-gluon contributions can effectively be dealt with. Obviously, this 
approach can be applied to D — > nK decay. 

The QCD factorization method can be applied to D — > -kit, nK and Kp decay, but we 
should calculate the contribution of power term 0( A< ^ D ). The power term includes the 
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contributions of soft-gluon effect, final state interaction and annihilation effects, since the 
power term in D — > nir, nK, Kp decay is larger than B — )■ nn, irK, Kp decay. We firstly 
considered D — > -kit decay in QCD factorization and light cone QCD sum rules method 
[7]. We found either the hard-gluon effect (0(a s ) correction) or the soft-gluon effect is 
small, and only found the calculation result of D° — > -n + it~ decay channel approaches 
the experiment data. It indicated that we should consider the contributions of final state 
interaction and annihilation effects in D — > tttt decay. In this paper, we apply the QCD 
factorization including light cone QCD sum rules method to study D — > irK, KK decay 
and obtain new results. In D° — > 7r°K° decay, we find both hard gluon and soft gluon 
contributions exceed the leading order largely, and the calculation result is accordance 
with the experiment data. In other decay channels, we should calculate all power terms, 
which include soft-gluon exchanges, final state interaction and annihilation effects, and 
then we can compare the calculation results with the experiment data. However, the final 
state interaction and annihilation effects haven't reliable method to calculate up to now. 
In our work, we calculate the leading order and a s corrections in QCD factorization, and 
the soft-gluon effects in the light cone QCD sum rules for the D — > irK, KK decay. In 
color-suppressed decay channels, we find the soft-gluon contributions are larger than the 
leading order contributions, and the calculation results are close to the experimental data 
for these decay channels. In color-allowed decay channels, the soft-gluon contributions 
are small and we should consider other power terms, such as the final state interaction 
and annihilation effects. 

2. D -> ttK^KK in QCD Factorization 

The low energy effective Hamiltonian for D° — > tt K° can be expressed as follows: 

tUff = ^OU^iMOiM + C 2 {p)0 2 {p)i (2) 

where Ci(p) are Wilson coefficients which have been evaluated to next-to-leading order. 
The four-quark operators Oi j2 

01 = (ud) v „ A (sc) V -A, 

2 = {u a d,3) v _ A (spc a ) v _ A , (3) 
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the Wilson coefficients evaluated at ti — m c scale are [8] 



d = 1.274, C 2 = -0.529. (4) 

In the following, we study D — > nK, KK decay with the QCD factorization approach. 
This method is similar to that for B —¥ -kit, ttK decay; see Ref. [5] for detail. As in Ref. 
[5], we obtain the QCD coefficients a$ at next-to-leading order (NLO) and 0(a s ) hard 
scattering corrections in naive dimension regularization (NDR) scheme. The coefficient 
a,i(7rK)(i = 1,2) are split into two terms: a,i(irK) = dij{nK) + a^n^K). They are given 
in Refs. [5, 14]. In D — > ttK,KK decay, the flavor structure is different from B decays. 
When we replace the index K and ir in the B decays' coefficients aiijiK){i = 1,2), we 
can obtain the coefficients a,i(nK)(i = 1,2) in D decays. They are 



n C 2/ C F a s C 2 C F na s 
a v = C X + - (1 + — K), aiJI = W ^-H„ K , (5) 



n , Cl n , C F a s C 1 C F na s 
a 2J = C 2 + Wc (l + —V K ), a 2JI = j^-^-H^. (6) 

Here N c = 3(/ = 4) is the number of colors (flavors), and C F — ^f^- is the factor of 
color. The functions in Eqs. (5) and (6) can be found in Ref. [5], which are 



V K = 12 In— - 18+ / g(xU K (x)dx, 
H Jo 

777 

V„ = 12 In— - 18+ / g(x)(f) n (x)dx, 
a Jo 



H 

1 — 2x 2 In a; 
g[x) = 3(- lax- in) + [2Li 2 (x) - (lnrr) 2 + (3 + 2ivr) hix-(x+±l- x)], 

U f°f R f 1 AC [ ldX A ( \ [ ld y\A f \, 2 ^KX 1 

h «k = 2 i7 D— v k f c\\ / —r~ d £ / —4>A X ) / —[My) + 



m 2 D F D ^ K (0) Jo £ Jo x ^ Jo y m c x ' 

H «* = mlF^(0) Jo —VJo Y Mx) Jo J [MV) + (7) 

where Li 2 (x) is the dilogarithm, fxifo) is the K (D meson) decay constant, m D is the D 
meson mass, F D ^ K (0) (F D_>7r (0)) is the D — >■ K (D — >■ ir) form factor at zero momentum 
transfer, and £ is the light-cone momentum fraction of the spectator in the D meson. 
H w k and Hk-k depend on the wave function 0£> through the integral Jo d£<Ao(£)/£ = 
bid/Ad = 6.23, with = (250 ± 75)MeV, \±k = 'rn 2 K /{m d + m s ), fi w = m%/(m u + m d ), 
m u = 3MeV, m d = 6MeV, m s = (150 ± 20)MeV, m c = l.SGeV, m n = 0.139GeV, 



m K = 0A94GeV. We take f n = (132 ± 0.26)MeV, f K = (170 ± 1.5)MeV, f D = (200 ± 
20)MeV, f Ds = (280 ± !8)MeV, F D ^(0) = (0.65 ± 0.10), F D ^ K (0) = (0.73 ± 0.07) [8], 
F D ^ K (0) = (0.82 ± 0.15) [15], a s (m c ) = 0.353, m D = 1.869GeV, m Ds = 1.968GeV, and 
the asymptotic wave functions (fix = (fin = 6x(l — x). 
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3. D — > ttK, KK in the light-cone QCD sum rules 



In the following, we calculate the soft-gluon contributions for D — > ttK, KK decay. 
Firstly, we calculate the soft -gluon effects of D° — > n !^ channel, and the calculation of 
other decay channels are similar as the channel. To estimate the soft-gluon corrections for 
D° — > 7r°~K° channel, it is useful to rewrite down the effective Hamiltonian with the help of 
the Fierz transformation. For example, applying the Fierz transformation to the operator 
2 = (uT^cj^T^d), we have the effective Hamiltonian relevant to the tree operators, 

Keff = ^v; s v ud [(cM + ^)om + 2CMOM], (8) 

where 

Oi = (uT^d)(sT^c), (9) 
In the above = 7/i (l - 75), Tr(\ a \ b ) = 25 ab , and 

2 = ^Oi + 20!. (10) 

First of all, we calculate the nonfactorizable matrix elements induced by the operator 
0\. As a starting object for the derivation of LCSR we choose the following vacuum-pion 
correlation function: 

F^\p,q,k) = -j d 4 xe-^J d%e^- k \0\T{j^(y)dmJ°(x)}\Aq))i (11) 

where jj^ ^ = 57^75^ and = m c ci'~f 5 u are the quark currents interpolating K° and D 
meson, respectively. The decomposition of the correlation function (11) in independent 
momenta is straightforward and contains four invariant amplitudes: 

Fj? 1} = (p - k) a F^ + qjf^ + k a F 2 (5l) + e^p^pf^. (12) 

In what follows only the amplitude F^°^ is relevant. The correlation function is calculated 
in QCD by expanding the T-product of three operators, two currents and 0±, near the 
light-cone x 2 ~ y 2 ~ (x — y) 2 ~ 0. For this expansion to be valid, the kinematical region 
should be chosen as: 

q 2 = p 2 = k 2 = 0, \(p - k) 2 \ ~ \(p - qf\ ~ \P 2 \ » A QCD 2 , (13) 
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where P = p — k — q. The correlation function (11) can be calculated employing the 
light-cone expansion of the quark propagator [6]: 



S(x,0) = -i{0\Tq(x)q(0)\0) 
T(d/2)x 

+ 16^2(^)^2-1 / M(i - v)xa^(vx) + va^(vx)x), (14) 

where = g s G a ^(\ a /2), which is the gluonic field strength and the soft gluon effects 
are from the term, d is the space-time dimension. Following the standard procedure for 
QCD sum rule calculation, we can obtain the hadronic matrix element of the operator 0\ 

A&\D° ->■ n°K°) = (^{-q^ip^O^ip - q)) 

= oj r 9 / dse^~ / ds 

7r 2 fKjDm D 2 Jo J ml 

e TP 1 Im sl Im s F^c D { s ^ s \ m2 D)^ ( 15 ) 



where Sq and Sq are effective threshold parameters for K and D meson. 

A straightforward calculation shows that only the twist-3 wave function (^(a,) and 
the twist-4 ones (p\\(ai), (p±(cxi), whose definitions can be found in Ref. [6], contribute to 
the invariant function F^ 0l \ The results are: 

F QCD — ^twZ + ftwi , UDJ 



with 



r (Oi) "T,c/3tt Z" 1 , , , , 



X 



(ml — (p — q) 2 (l - a 1 ))(-P 2 va 3 — (p — k) 2 (l - va 3 )) 
x [(2 - v)(q ■ k) + 2(1 - v)q ■ (p - k)](q ■ [p - k)). (17) 

We get the same calculation results as Ref. [6] for Eq. (17) when it be substituted for 
m c — > m b , but the twist-4 contribution hasn't been showed in Ref. [6]. Now, we give the 
invariant amplitude from the twist-4 term. 

<o.) _ KU f 1 , f 1 (4tt-6)(p-fc)g 



m c — (p ~ Q + Q a i) 2 (p — k — qvaz) 2 



m 2 f n f 1 f 1 (2pq — 2vqk)(p — k)q 



dv J daida3$i(ai : a?) 



2ir 2 Jo J ' \m 2 — (p — q + qa>i) 2 ] 2 (p — k — qva^) 2 

J dv J da 3 § 2 (a 3 ) 



m 2 c U f 1 dy f da ^j a \ I (2pq - 2vqk)(p - k)q 

2ir 2 Jo J [m 2 — (p — qa^) 2 ] 2 (p — k — qva^) 2 

4 ™> 2 cfv f 1 j.. ..2 [ j„. ^ („. \ ' 1 ' [( P -k)qf 



■ J dv2v 2 J da 3 § 2 (a 3 ) 



2ir 2 Jo J pq[m 2 c — (p — gc^) 2 ] (p — k — qva 3 ) 4 



, [ 1 dv(2v-2)v [ da 3 $ 2 (a 3 )^-± ^ [(P (18) 

27r 2 Jo 7 m 2 — (p — qa 3 ) 2 (p — k — qva 3 ) 4 

In Eq. (17) and (18), we make use of the following nonlocal operator matrix elements: 

(0\u(0)a^ l5 G a p(vy)u(x)\7r°(q)) 

= i^[(q a q»gpis - w^c.v) - (q a qug^ - qpqvg^)} 

x J Dam^^e-^+y^ , (19) 
(0\u(0)i llx G a p(vy)u(x)\ir o (q)) 

- Da^^,^-^^, (20) 



with Dai = daida 2 da 3 S(l — ot\ — a 2 — 013). Finally, the LCSR for the D° — > n^TT matrix 
element of the operator 0\ from the soft-gluon exchange is obtained by applying to the 
duality approximation and Borel transformation in the D channel. The result can be 
written as: 

= im n\ — 7tt- / dse m 1 )( — j- I — c M um> 2 

rrnjto [»dv n dv \a~ n \ 

X / ^3ir(l -U,U-V,V) + / 7r / 3<^_l(1 -U,U- V,V) 

u Jo v Jo V 

uM' 2 > u l + M uM' 2 ] u 2 lh [Zi) 
where Uq = m 2 / 1 ', and the following definitions are introduced: 

<9$i(-u;,f) 



dw 
d$ 2 (v 



(f±(w, 1 — w — v, v) + y?||(w, 1 — w — v , i>), 



(> =$i(l -«,«). (22) 
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The asymptotic forms of the pion distribution amplitudes in Eqs. (21)- (22) are given by 
[6]: 

y?37r(tti) = <P3k(<Xi) = 360aia2«3, 

<P±(<Xi) = I05 2 al(l - a 3 ), 

(f>\\(ai) = -A05 2 a 1 a 2 a 3 . (23) 

We find our calculation result (Eq. (21)) has a little difference with Eq. (30) in Ref. 
[6]. In Eq. (21), the final term includes function (^k — 2)^|^-, and this function is 
corresponding to (^gk — jjn ~ 1) $2 i"^ in Ref. [6] , which includes variation s. We think the 
function including variation s should appear in the first integral /q° dse^ 1 , and it should 
not be in the second integral • • • du, so that the matrix element A is independent of 
variation s. Our calculation result A\ (Eq. (21)) is a constant. In Ref. [6], the matrix 
element A (Eq. (30)) is a function of variation s. It should print error in Ref. [6]. 

We write the decay amplitudes of D° — > 7T°K°, which include factorization and non- 



factorization parts: 



M f+as (D» -> tt u 10 = i^V; s V ud f K F D ^(0)(m 2 D - m*)a 2 , (24) 



M nfg (D° -+ n°K°) = ^PlG F V: s V ud C x A x . (25) 

The matrix elements of other decay channels, can be calculated as similar as D° — >■ 
7r°7C° channel, and we can write down them directly as follows: 

A 2 = A^\D° ^n + K-) 

DK 4n 2 fJo JK 2f D m%Ju» u 

x / — <P3k{1 - u,u- v,v) + } K — 3</?_l(1 -u,u- v,v) 

u Jo V Jo V 

-l)*!fi^0] +/jr (^ I -2)^M]), (26) 



uM' 2 ' u J J uM' 2 ' u 2 



. 1 rs£ rtzl m 2 c /-i d ' 

DA A^f K Jo A 2f Ds m% 3 Ju» u 

vI m c /3X f u dv ( udv \o~ a \ 

x / — <P3k{i ~ u,u- v,v) + f K — 3</?_l(1 - u,u- v,v) 

u Jo v Jo V 



ml $i(l — u,v) ml $ 2 (" 



!) " - ' - 1 + /tt - 2)^]), (27) 



'tiM' 2 ; u J J uM' 2 ' u 2 
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and 



m 2 f . f 1 du ^Ds m? 

= ilTLn ( I dSe«^)( — j— I CM' 2 uM' 2 



A A = A( 02 \D+ — >■ nK) 

^<m c f 3K [ u dv ( udv \o~ a \ 

x / — (P3k{1 -u,u-v,v) + J k / — 3</?_l(1 -u,u-v,v) 

u Jo v Jo v 

K_ *i(l^Mi ]+M < _ Mf)]), (28) 
mm' u um' m z 

where 2 = (wr M ^c)(sr M 4^-<f) and is effective threshold parameters for 7r meson. 

In the following, we give the decay amplitudes of other channels, which include fac- 
torization and non-factorization parts: 

M f+as (D° n+K-) = i 9£ V ; s V ud UF D - K (0)(ml - m 2 K )a u (29) 

M nfg (D° -> 7T + K~) = VtGpV'Vrt^. (30) 

M /+tt .(D + ->7r + 7r°) = i^^Kd^F^W^-m^a! 

+/^ D -(0)K- m ^a 2 ], (31) 



M n/g (L>+ tt+TT) = v^Gf^KdCdAx + C 2 A 2 ). (32) 
M f+as (Di ^ + ^°) = ^K;K^F M (0)(m^ - m^)a 2 , (33) 

^n/ ff (£> s + ^ + ^°) = "V^GF^KdCiAg. (34) 
M f+as (Dt tt+AT ) = »|w/ fl '^(0)K - <K (35) 

M n/ff (L>+ Tr+iT ) = y/2G F V*y vd C 2 A 4 . (36) 

where the amplitude Mf +CCs = Mf + M Qs represents the sum of the leading order fac- 
torization Mf and non-factorization M as from the hard-gluon exchanges, the amplitude 
M n f g is the non-factorization parts from soft-gluon exchanges. The total amplitude M is 
the sum of Mf +as and M n f g . 
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4. Numerical calculation 



In the numerical calculations we take Sq = 0.7Gey 2 [6], sff = 1.2GeV 2 [9] and = 



(6 ± l)GeV 2 [lO]. In fx c = ^m 2 D -m 2 c « l.ZGeV , f 37V (fi c ) = 0.0035Ge\/ 2 , <5 2 (/i c ) = 
0.19GeV 2 are nonperturbative parameters in light-cone wave functions [10], and the CKM 
matrix are V* d = V cs = 0.9734 + 0.9749 and V cd = 0.227 [11]. Having fixed the input 
parameters, one must find the range of the values M' 2 and M 2 for which the sum rules (Eq. 
(21)) is reliable. At the interval M' 2 = 8-12GeV 2 and M 2 = 6-15GeV 2 , we find the value 
of Ai (Eq. (21)) is quite stable. The D meson life time r(D°) = (4.12 ± 0.027) x 10~ 13 s, 
T (D + ) = (1.05 ±0.013) x 10- 12 s, t(D+) = (4.9 ±0.09) x 10~ 13 s. In the D rest frame, the 
two body decay width is 

T(D -> P 1 P 2 ) = ^-\M(D -> P^f^l, (37) 

where P\ and P 2 are two pseudoscalar meson (tt and K), and the momentum of the Pi 
meson is given by 

I pi [( m2 D ~ ( m Pi + m P 2 ) 2 )K ~ ( m Pi - m P2?)]^ /ooA 

|P| - 2^ ' (38) 

The corresponding branching ratio is given by 

Br(D -+ Pl P 2 ) = ^V^ - (39) 

J- total 

where T tota i denotes the total decay width of D meson. The total decay width of one 
meson is related to its mean life time r by T tota i = Ti/t. With the above parameters 
and formulae, we can get the branching ratios of D — > irK, KK decay obtained in some 
approaches with that of the experiment. 

Table 1: The branching ratios of D — > nK, KK decay obtained in some approaches 
together with experimental result. 



Decay channel 


NF 


QCDF 


QCDF+SGE 


Experiment 


D° -> TT ^ 


2.4 x 10~ 3 


(3.66 ±0.55) x 10~ 2 


(2.20 ±0.11) x 10~ 2 


(2.28 ±0.22) x 10~ 2 


D {> ->■ TT+K' 


5.63 x 10~ 2 


(7.18 ± 1.01) x 10~ 2 


(6.15 ±0.40) x 10~ 2 


(3.80 ±0.09) x 10~ 2 


D+ ->■ TT+Tf 


6.55 x 10~ 2 


(1.66 ±0.23) x 10~ 2 


(2.70 ±0.17) x 10~ 2 


(2.77 ±0.18) x 10~ 2 


D+ ->■ K + K u 


2.15 x 10~ 3 


(4.27 ±0.68) x 10~ 2 


(3.86 ±0.26) x 10~ 2 


(3.6 ± 1.1) x 10~ 2 


D+ ->■ n+K° 


3.49 x 10~ 3 


(5.64 ± 0.97) x 10~ 3 


(5.07 ±0.35) x 10~ 3 


< 8 x 10~ 3 
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The branching ratios of D — > ttK, KK decay channels are presented in Table 1, where 
the second column is the result of naive factorization (NF) and the total amplitudes 
Mf +(Xg in Eqs. (24), (29), (31), (33) and (35) corresponding to different decay channels 
of D — > irK, KK, a± and a 2 are calculated in the leading order, i.e. the parameters 
o-i — Ci + ^ and a 2 = C 2 + the third column is the result of QCD factorization 
(QCDF), the amplitudes are also calculated by M f+OLs in Eqs. (24), (29), (31), (33) and 
(35) but ai and a 2 are calculated by QCD factorization approach from Eqs. (5)-(7), which 
include the leading order and 0(a s ) corrections, the fourth column is the result of QCD 
factorization including soft-gluon exchanges (QCDF+SGE) which are our results, the 
total amplitude is the sum of Mf +olg and M ng f in Eqs. (24)-(25) and (29)-(36 ), and the 
final column is the experimental data [12]. From Table 1, we can find that the prediction 
of naive factorization is far from the experimental data and the QCD factorization method 
have improved the calculation results. In our approach (QCD factorization including soft- 
gluon effects), the calculation results are close to the experimental data in D° — > tt K°, 
D + — > 7i + K° and D+ — > K + K° decay channels. We find the soft-gluon corrections 
are rather large, which are larger than the leading order contributions in D° — > 7r°X° 
and D+ — > K + TC° decay channels. For D° — > n + K~ decay, the results from the three 
approaches do not agree with the experimental data, and we think the reason is that we 
only calculate the soft-gluon contributions and do not consider the final state interaction 
and the annihilation effects in the power term 0(AQ C D/ m c)- The theoretical uncertainties 
in the table 1 are estimated by the some parameters. We can give the uncertainties 
from QCD factorization and soft-gluon effects, respectively. In QCD factorization, the 
important uncertainties are from the input parameters: form factors F D ~* K , F Ds ^ K , 
pD^-K ^ decay constants f D) f Ds) f K) fn and the parameters ai and a 2 - From Eqs. (5)- 
(7), we can find the uncertainties of ai and a 2 are from: (1) The Wilson coefficients: 
The coefficients C\ and C 2 are uncertain, and they are in the ranges of: 1.216 ~ 1.274 
and —0.415 ~ —0.53 respectively. (2) The vertex corrections Vm{M = tt,K): we need 
to input the light-cone distribution functions <Pm{M = it,K). We take the asymptotic 
form of pion and kaon light-cone distribution functions <P-k{x) = <Pk{x) = Qx(l — x), 
but the accurate form should be Gegenbauer polynomials. (3) The hard-scattering terms 
H K tt and H^k'- we need to calculate the moment j-^oiO = the value of Ad at 
present is uncertain, a typical range being Ad = (250 ± 75)MeV. The uncertainties from 
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QCD factorization less than 18%. In soft-gluon effects, the uncertainties are from the 
input parameters: the decay constants fr>, fu a , fx, f n , the Wilson coefficients C\ and 
C 2 , and the Borel parameters M' 2 and M 2 . The uncertainties from soft-gluon effects 
less than 7%, and the total uncertainties less than 25%. In D — y ttK, KK, we find 
the QCD factorization (QCDF) and QCD factorization including soft-gluon exchanges 
(QCDF+SGE) have improvement on the calculation results. So the QCD factorization 
method can be applied to D — y nK, KK decay, but the power term 0(Ag C . D /m c ) should 
be included. In order to get accurate calculation results in D — > ttK, KK decay, the soft- 
gluon, the final state interaction and the annihilation effects should be calculated, and the 
calculation results can be improved. Now, the final state interaction and the annihilation 
effects have many models, but there is not a reliable method. So, it is necessary to study 
the D — y nK, KK decay further. 

5. Summary 

In D — y nK, KK decay, We find the prediction of naive factorization (NF) is far from 
the experimental data, the QCD factorization (QCDF) and QCD factorization including 
soft-gluon exchanges (QCDF+SGE) have improvement on the calculation results. In 
color-suppressed decay channels, such as D° — y n°K° and Df —y K + K° decay, the soft- 
gluon corrections are rather large, which are larger than the leading order contributions. 
In color-allowed decay channel D + — y n + K°, the soft-gluon corrections are rather large 
also, but the soft-gluon corrections are small in the color-allowed decay channel D° —y 
n + K~ . From Refs.[6, 13], we can find the soft-gluon contributions are different in different 
B — y 7T7T decay channels. For example, in B° — y tt + tt~ decay, the soft-gluon effects are 
smaller than the hard-gluon contributions because of the small Wilson coefficient C 2 in the 
decay amplitude. However, the soft gluon effects and the hard gluon are on the same order 
in B° —y tt°-k° decay since the Wilson coefficient C\ is large in this decay amplitude. In 
D —¥ irK, KK decay, we can obtain the similar results. The QCD factorization method 
can be applied to D — y irK,KK decay, but the power term 0(Ag CD /m c ) should be 
included. In order to get accurate calculation results in D — y irK, KK decay, the soft- 
gluon, the final state interaction and the annihilation effects should be calculated, and 
then the calculation results can be improved. 
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